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Effective low-energy Hamiltonians for several different families of iron-based su- 
perconductors are compared after deriving them from the downfolding scheme 
based on first-principles calculations. Systematic dependences of the derived model 
parameters on the families are elucidated, many of which are understood from the 
systematic variation of the covalency between Fe-3d and pnictogen-/chalcogen-p 
orbitals. First, LaFePO, LaFeAsO (1111), BaFe 2 As 2 (122), LiFeAs (111), FeSe and FeTe (11) 
have overall similar band structures near the Fermi level, where the total widths of ten-fold 
Fe-3d bands are mostly around 4.5 eV. However, the derived effective models of the ten-fold 
iron-3d bands (d model) for FeSe and FeTe have substantially larger effective onsite Coulomb 
interactions U ~ 4.2 and ~ 3.4 eV, respectively, after the screening by electrons on other bands 
and after averaging over orbitals, as compared to ~ 2.5 eV for LaFeAsO. The difference is 
similar in the effective models containing p orbitals of As, Se or Te (dp or dpp model), where U 
ranges from ~ 4 eV for the 1111 family to ~ 7 eV for the 11 family. The exchange interaction 
J has a similar tendency. The family dependence of models indicates a wide variation rang- 
ing from weak correlation regime (LaFePO) to substantially strong correlation regime (FeSe). 
The origin of the larger effective interaction in the 11 family is ascribed to smaller spread of 
the Wannier orbitals generating larger bare interaction, and to fewer screening channels by 
the other bands. This variation is primarily derived from the distance h between the pnicto- 
gen/chalcogen position and the Fe layer: The longer h for the 11 family generates more ionic 
character of the bonding between iron and anion atoms, while the shorter h for the 1111 family 
leads to more covalent-bonding character, the larger spread of the Wannier orbitals, and more 
efficient screening by the anion p orbitals. The screened interaction of the d model is strongly 
orbital dependent, which is also understood from the Wannier spread. The dp and dpp models 
show much weaker orbital dependence. The larger h also explains why the ten-fold 3d bands 
for the 11 family are more entangled with the smearing of the "pseudogap" structure above the 
Fermi level seen in the 1111 family. While the family-dependent semimetallic splitting of the 
bands primarily consists of d yz /d zx and d x 2_ y 2 orbitals, the size of the pseudogap structure 
is controlled by the hybridization between these orbitals and d xy /d 3z 2_ r 2: A large hybridiza- 
tion in the 1111 family generates a large "band- insulating" -like pseudogap (hybridization gap), 
whereas a large h in the 11 family weakens them, resulting in a "half-filled" like bands of or- 
bitals. This may enhance strong correlation effects in analogy with Mott physics and causes the 
orbital selective crossover in the three orbitals. On the other hand, the geometrical frustration 
t! jl, inferred from the ratio of the next-nearest transfer t' to the nearest one t of the d model 
is relatively larger for the 1111 family than the 11 one. The models comprehensively derived 
here may serve as a firm starting basis of understanding both common and diverse properties 
of the iron-based superconductors including magnetism and superconductivity. 

KEYWORDS: first-principles calculation, effective Hamiltonian, downfolding, constrained RPA method, 
LaFePO, LaFeAsO, BaFe2As2, LiFeAs, FeSe, FeTe, oxypnictide, oxychalcogenide, high- 
temperature superconductivity 



1. Introduction 

Recently a new class of superconducting compounds 
including iron element has been discovered. 1 ) In all the 
cases of this class, Fe-3d conduction electrons are likely 
to form Cooper pairs and responsible for the supercon- 
ductivity. However, the mechanism of superconductiv- 
ity is not well understood and is under extensive de- 
bates. In the family with ZrCuSiAs-type structure (called 
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1111 hereafter), SmFeAs(0,F) has shown the highest su- 
perconducting critical temperature T c ~ 56 K (ref. 2) 
when fluorine is substituted for ~ 20% of oxygen as elec- 
tron doping, while BaFe2As2 with ThCr2Si2-type struc- 
ture (called 122) has indicated the highest T c ~ 38 
K, when potassium is substituted for ~ 40% of Ba as 
hole doping. 3 ^ There exist another simpler compounds 
LiFeAs and NaFeAs (called 111) reported as the PbFCl- 
type tetragonal structure, indicating T c ~ 18 K. 4 ~ 6 ' An- 
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other family of binary compounds FeSe^Tei-^ (called 
11) also shows superconductivity at T c higher than 10 
K (refs. 7, 8) and has reached 37 K under pressure (7 
GPa). 9 ' 

It is highly desired to understand what are common 
and what are family dependent, from detailed electronic 
structure of these four families, for the purpose of es- 
tablishing the basis for revealing mechanisms of super- 
conductivity and magnetism. In this report, we present 
effective low-energy models of these families derived from 
first principles and compare them. 

A common known feature of iron-based superconduc- 
tors is that the antiferromagnetic order appears close 
to the superconducting region except for the 111 fam- 
ily. The mother compound of the 1111-type, for exam- 
ple, LaFeAsO shows antiferromagnetic long-range order 
of the stripe type below Tjv ~ 130 K with the Bragg 
point at (71", 0) in the extended Brillouin zone (corre- 
sponding to (7r, 7r) in the reduced Brillouin zone). 10 ) The 
antiferromagnetic ordered moment ~0.36 hb as com- 
pared to the nominal saturation moment 4 \ib for the 
high-spin 3c? 6 state is unexpectedly small, implying large 
quantum fluctuations arising from electron-correlation 
effects or dominating itinerancy with subtlety of com- 
peting ground states. On the other hand, the 122-type 
(BaFc2As2) shows a relatively large ordered moment 
~ 1.1 [i b (refs. 11 and 12) and the 11-type (FeTe) in- 
dicates even larger ordered moment ~ 2.25 \ib at a dif- 
ferent Bragg point, (tt/2,tt/2) in the extended Brillouin 
zone. 8 ' Even for the 1111 family, NdFcAsO shows larger 
ordered moment (~ 0.9 /is), although the moment is 
apparently reinforced by Nd moment (~ 1.55 /ib)- 13 ^ 

Conventional density-functional calculations with the 
local density approximation (LDA) or the generalized 
gradient approximation (GGA) have clarified that bands 
originating from ten-fold degenerate iron-3c£ orbitals in 
a unit cell containing two Fe atoms are close to the 
Fermi level. The LDA calculations of the 1111-type, 14 ~ 20 ' 
122-type, 21 ' 22 ' 1 11-type, 23 ' and 11-type compounds 24 ' 25 ' 
show a very similar band structure for all of the above 
families, where small electron pockets around M point 
and hole pockets around Y point lead to a semimetal- 
lic Fermi surfaces. The local spin density approximation 
(LSDA) also commonly predicted the antiferromagnetic 
order for mother materials. 16,17,19 ' The stripe-type anti- 
ferromagnetic order is correctly reproduced for the 1111- 
type. 17 ' 19 - 1 However, the calculated ordered moment ob- 
tained so far is large and ranges between 1.2 and 2.6 
Mb, 16 ' 17 ' 19 ' 26 ' 1 in contrast to much smaller ordered mo- 
ment discussed above. It is unusual to observe the or- 
dered moment smaller than the LSDA result. On the 
other hand, the bicolincar order for FeTe is reproduced in 
the LSDA with more or less consistent ordered moment 
(~ 2.25 hb) with the experimental results. 25 ' Broad peak 
structures of magnetic Lindhard function calculated by 
using the LDA/GGA Fermi surface suggest competitions 
of several different ordering tendencies. 18 ' 27 ~ 30 ' 

The role of electron correlation based on the realistic 
grounds is under a strong debate. 31 ~ 34 ' Relatively small 
fraction of the Drude weight 35-38 - 1 together with bad 
metallic (or semiconducting) behavior 1 ' 39 ' indicate sub- 



stantial electron correlation effects. Antiferromagnetic 
orders and fluctuations themselves revealed by the nu- 
clear magnetic resonance and other probes near the su- 
perconducting phases also indicate, in any case, some 
electron correlation effects play a role. 40 ' 41 ' Diversity of 
the ordered moment ranging from 0.36 fiB to 2.26 /ib is 
remarkable in terms of the similar band structure with 
scmimctallic small pockets of the Fermi surface. Ordering 
vector of the antiferromagnetic order introduced above 
depending on the compounds further suggests that the 
correlation effects and its subtlety arc beyond the sim- 
ple nesting and weak coupling picture. In fact, recent 
fluctuation exchange calculation suggests that the self- 
energy effect with subtle multiband structure near the 
Fermi level cast a serious suspicion on the validity of the 
nesting picture. 42 ' 

Angle resolved photoemission spectroscopy 43,44 ' has 
shown some correspondence to the LDA result of Singh 
et aZ. 21 ' Ye-2p core-level spectra of X-ray photoemission 
suggest rather itinerant character. 45,46 ' On the other 
hand, some role of moderate electron correlations has 
also been claimed. 47 ' 48 ' For FeSe, soft-Xray photoemis- 
sion results 49,50 ' appears to show a deviation from the 
LDA results suggesting a splitting of the coherent band 
near the Fermi level from the incoherent part arising from 
the correlation effect as we discuss in Sec. 4. 

In the superconducting phase, even the pairing sym- 
metry itself is highly controversial and no consensus 
has been reached. Although nodeless superconductivity 
is suggested, 51 ~ 53 ' temperature dependence of nuclear- 
magnctic-relaxation time T\ below T c roughly scaled by 
T~ 3 without the Hcbcl-Slichtcr peak implies unconven- 
tional superconductivity driven by nontrivial electron- 
correlation effects. 40 ' For example, orbital dependent 
gaps with sign-changing and fully-gapped s± symmetry 
has been proposed. 41 ' The gradual suppression of the su- 
perconducting transition temperature by Co doping into 
the Fe site was reported to be explained by strong antifer- 
romagnetic fluctuations near the metal-insulator bound- 
ary producing an effect on the ,s-wavc singlet pairing 
without the sign change. 54 ~ 56 ' Although overall experi- 
mental results suggest noticeable correlation effects, real- 
istic roles of electron correlations on theoretical grounds 
are not well established and controversial. 

The effective electron Coulomb repulsion of the 1111 
family for the models of the ten bands of Ye-3d orbitals 
has been estimated by Nakamura et al. from first prin- 
ciples by applying the downfolding scheme to eliminate 
other band degrees of freedom 32 ' along the line of the 
three-stage scheme. 57 ~ 59 ' The ratio of the Hubbard on- 
site interaction U ~ 3 eV to the nearest neighbor transfer 
t ~ 0.3 cV in the downfoldcd model is estimated to be 
U/t ~ 10 with the fivefold orbital degeneracy, indicating 
a moderately strong correlation. This identification has 
also been supported for the case of the 122-type. 60 ' Here 
it should be noted that the effective Coulomb interaction 
estimated in these works as well as in the present paper 
for the low-energy model is not directly the same as the 
interaction derived in experimental probes such as the 
X-ray photoemission. 45,46 ' This is because, in the model 
parameters here, the screening arising from the polar- 
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ization within the low-energy degrees of freedom (Fe-3d 
bands in ref. 32) is excluded as we describe in the next 
section, because this screening effect should be consid- 
ered when the low-energy model is solved. On the other 
hand, in the experiment, the interaction effect appears 
as a whole consequence after the full screening. We will 
further discuss this issue in Sec. 4. 

When we consider the puzzling diversity of magnetic 
and transport properties among the 1111, 122, 111 and 
11-types in spite of the apparent similarity of the band 
structure by the LDA, it is crucially important to eluci- 
date the origin of the difference and diversity from a uni- 
fied first-principles calculations by taking into account 
electron correlations properly. In this paper, we extend 
the work for the 1111 family 32, 61 > 62 ) and derive effective 
low-energy models of the 1111 as well as 122, 111 and 
11 families on a unified grounds. We further classify and 
compare effective models toward the comprehensive un- 
derstanding of the electron correlation effects. This com- 
parison on the diversity is also important for the under- 
standing of the superconducting mechanism. A key quan- 
tity is the pnictogen/chalcogen height, h. It was pointed 
out in the early stage that the electronic band struc- 
ture is altered significantly by changing /i. 15,63 ) There 
is also a previous work claiming that the spin and the 
charge susceptibility are sensitive to /i. 30 ) The present 
study reveals that the strength of electron correlation is 
determined by the spatial extent of the Wannier orbitals 
and the strength of screening effect, and both of them 
are affected by h. 

In Sec. 2 we describe our method. Sec. 3 describes 
the derived effective models for LaFcPO, LaFeAsO, 
BaFe2As2, LiFeAs, FeSe and FeTe. We present effective 
models both for the ten-band model for the iron-3e? Wan- 
nier orbitals (d model) and the model including p orbitals 
of P, As, Se or Te (dp or dpp model). Sec. 4 is devoted 
to summary and discussions. 

2. Method 

We derive the low-energy models by a combined con- 
strained random-phasc-approximation (cRPA) (ref. 57) 
and maximally localized Wannier function (MLWF) 
(refs. 64 and 65) method. This combination has been 
recently developed and successfully applied to the 3d 
transition metals, 66,67 ) their compounds 66 -* including 
LaFeAsO, 32,61,62 ) organic conductors, 68 ) and zeolites. 69 ) 

The first step of the method is a standard band struc- 
ture calculation in the framework of density functional 
theory (DFT). 70, 71 ) We then choose target bands around 
the Fermi level and extract them following the MLWF 
procedure, which defines the one-body part of the low- 
energy model. The transfer integral is obtained by tak- 
ing the matrix element of the Kohn-Sham Hamiltonian, 
Hks, in the MLWF basis, 



tmn (-^) 



(1) 



where 0„jj('") is the MLWF centered at the site R for 
the n-th orbital. 

To evaluate the effective interaction parameters, par- 
tially screened Coulomb interaction at zero frequency, 
W r (r, r'; u> = 0), is calculated in the cRPA with the con- 



straint that screening channels inside the target bands 
are cut out. This constraint is imposed to avoid dou- 
ble counting of the screening effects; they are consid- 
ered later when the derived model is solved. The effective 
Coulomb interaction U and exchange interaction J are 
orbital dependent. Their matrix elements are given by 

U mn (R) = (0mO'/ l mo|W / r|0„_R0 n _R) , (2) 

J mn (R) = ((f) m o(t)no\W r \(j) nR (j) mR ) , (3) 
(4>i4>j\Wr\h4>i) = J [ 4>*(r)4>j(r)W r (r,r';u; = 0) 

x<f)* k (r')<f)i(r')drdr' . (4) 

One problem in the cRPA method is that, in the case 
of entangled band structure, it is not clear which screen- 
ing process is to be excluded. This is indeed the case in 
all the materials studied in this work except for FeSe. Ex- 
tension of the cRPA technique for entangled band struc- 
ture has been proposed very recently. 67 ) In this method, 
the Hilbert space is divided into two parts; low-energy 
space spanned by the MLWFs and the rest of the space. 
Neglecting hybridization between the two spaces and di- 
agonalizing the Kohn-Sham Hamiltonian in each space, 
disentangled band structure is obtained. Screening chan- 
nel inside the target space is well-defined in the disen- 
tangled band structure, hence the cRPA calculation can 
be done without ambiguity. The present study serves as 
a good application of this technique. 

The band structure calculation is based on the 
full-potential LMTO implementation. 72 ) The exchange- 
correlation functional is the local density approximation 
of the Cepeley- Alder type 73 - 1 and spin-polarization is ne- 
glected. The cRPA calculation uses a mixed basis con- 
sisting of products of two atomic orbitals and interstitial 
plane waves. 74 ) The self-consistent LDA calculation is 
done for the 12 x 12 x 6 fc-mesh, and 20 x 20 x 10 k 
points are sampled for the density of states. Both the 
partially screened Coulomb interaction and the MLWF 
setup use the 4x4x4 mesh. More technical details are 
found elsewhere. 66 -* 

The cRPA calculations were also performed with an- 
other ab initio band-structure code based on plane-wave 
basis set, Tokyo Ab initio Program Package, 75 * 1 for crit- 
ical comparisons with the FP-LMTO results. Density- 
functional calculations with LDA within the parame- 
terization of Perdew-Wang 76 ) were performed with the 
Troullier-Martins norm-conserving pseudopotentials 77 ) 
in the Klcinman-Bylander representation. 78 ' Iron pseu- 
dopotential was constructed under the reference config- 
uration (3d) 7 (4s)°- 8 (4p) - by employing the cutoff ra- 
dius for the 3d state at 1.3 Bohr and for 4s and Ap states 
at 2.1 Bohr, with supplementation by the partial core 
correction of cutoff radii of 0.6 Bohr. The cutoff energies 
in wavefunctions and charge densities were set to 100 Ry 
and 900 Ry, respectively, and a 5 x 5 x 5 fc-point sampling 
was employed. The polarization function was expanded 
in plane waves with an energy cutoff of 20 Ry and the 
total number of bands considered in the polarization cal- 
culation was set to 130. The Brillouin-zone integral on 
wavevector was evaluated by the generalized tetrahedron 
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method. 79 -* The additional terms in the long-wavelength 
polarization function due to nonlocal terms in the pseu- 
dopotentials were explicitly considered following ref. 80. 
A problem due to the singularity in the Coulomb inter- 
action, in the evaluation of the Wannier matrix elements, 
U mn (R) and J mn (R) in Eqs. (2) and (3), was treated in 
the manner described in ref. 80. We checked that these 
conditions give well converged results. 

A dependence of the resulting screened onsite parame- 
ters on cutoff radius of the iron pscudopotcntial was care- 
fully checked through calculations with different choices 
of the cutoff radii of 1.3, 1.7, and 2.1 Bohr. We found that 
these different choices make no discernible difference in 
the resulting values and the difference is less than 0.1 eV 
at maximum; for example, for U xy of LaFeAsO, the value 
is 3.14 eV for r c = 1.3 Bohr, 3.20 eV for r c = 1.7 Bohr, 
and 3.17 eV for r c = 2.1 Bohr. Notice that differences in 
the values are not necessarily monotonic. 

Now the polarization effects from the other bands far 
from the Fermi level are considered in the cRPA yielding 
the screened Coulomb interaction of the target bands. 
The degrees of freedom for the other bands are elimi- 
nated, leaving the low-energy degrees of freedom within 
the target bands only. Thus obtained low-energy effective 
Hamiltonian has the form 

K = y^y^ ^ t mn (Ri - Rj)a1 n a a jm 

(7 ij nm 

+— / ] / ] s U mn (Ri — Rj)a in a P l n a P m a in 

ap ij nm 

~\~Jmn{.Ri — Rj ) (^in a jm a in a jm^ a in a in a %i a jm^ i (5) 

where (af n ) is a creation (annihilation) operator of 
an electron with spin a in the nth MLWF centered on 
Fe atom at R4. In Sec. 3 the target bands left in the low- 
energy bands are either Fe-3d bands (d model) or Fe-3d 
as well as p bands of pnictogen and oxygen (or chalcogen) 
(dp/dpp model). 

3. Results 

3. 1 Band structure and density of states 

Figures 1 and 2 show the band structures and densities 
of states, respectively, of LaFcPO, LaFeAsO, BaFc2As2, 
LiFcAs, FeSe and FeTe in the experimental geometry 
(Table I). The overall feature of the band structure 
is common in all the compounds, reflecting the com- 
mon existence of the Fe layer sandwiched by the pnicto- 
gen/chalcogcn atoms as is illustrated in Fig. 3. In Fig. 1 
we see entangled ten bands having strong Fe-3<i character 
located near the Fermi level. 14 ~ 25 ) The bandwidth ranges 
from 4.4-4.6 eV (LaFeAsO, BaFe 2 As 2 , FeSe, FeTe) to 4.9- 
5.1 eV (LiFeAs and LaFePO). The systems are metallic 
in the LDA with electron pockets around the M-A line 
and hole pockets formed around the T-Z line. Below the 
d band are three states, which are mainly of P-/As-/Se- 
/Tc-p character. 

Looking at the density of states, we find that the d 
band in the 1111 family has a dip (pscudogap) at the 
energy roughly 0.5 eV higher than the Fermi level. In 
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Fig. 2. Density of states of six different iron-based superconduc- 
tors obtained by DFT-LDA. Number of states is counted for 
one- half formula unit in BaFe2As2, and for one formula unit in 
the other materials. Energy is measured from the Fermi level. 



Table I. Lattice parameters used in the present work. Here, h is 
the distance between the pnictogen/chalcogen atom and the Fe 
plane as is illustrated in Fig. 3. 





a (A) 


c(A) 


h(A) 






LaFcPO 


3.9636 


8.5122 


1.1398 


ref. 


81 


LaFeAsO 


4.0353 


8.7409 


1.3216 


ref. 


1 


BaFe2As2 


3.9625 


13.0168 


1.3602 


ref. 


3 


LiFeAs 


3.7764 


6.3568 


1.5075 


ref. 


5 


FeSe 


3.7738 


5.5248 


1.4652 


ref. 


82 


FeTe 


3.8123 


6.2517 


1.7686 


ref. 


8 



other families, the pseudogap is not clear, and it is com- 
pletely smeared out in FeTe. The partial density of states 
and the occupation number resolved by the MLWF in 
the d model are shown in Fig. 4 and Table II, respec- 
tively, for LaFeAsO, FeSe, and FeTe. The occupation 
number is the largest for the d^ z 2_ r 2 orbital which is 
nearly | -filling in all the compounds. Hereafter we ab- 
breviate the 3d orbitals such as d 3z 2_ r 2 as 3z 2 — r 2 , un- 
less confusions occur. The x 2 — y 2 and yz/zx orbitals 
are roughly half-filling, although the weight of the for- 
mer increases as we move from LaFeAsO to FeSe and 
FeTe. It turns out that the x 2 — y 2 and yz/zx orbitals 
are the primary origin of the evolution of this pscudo- 
gap structure. The pDOS for these orbitals are strongly 
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Fig. 1. Electronic band structures of six iron-based superconductors obtained by DFT-LDA. The K\ — K§ points in BaFe2As2 are 
Ki = ^(|,0,0),K 2 = &.(L,l,0),K 3 = ^(0,0,^),iC 4 = ^(|,0,^),^5 = respectively. Energy is measured from 

the Fermi level. 
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Fig. 4. Color: (Left panel) Partial density of states of LaFeAsO, FeSe and FeTe resolved by the MLWF in the d model. Energy is 
measured from the Fermi level. Pseudogaps seen in yz/zx and x 2 — y 2 bands for LaFeAsO around E ~ 0-0.5 eV are filled in FeTe. 
(Right panel) Partial density of states without interorbital hybridization is plotted for comparison. 
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Fig. 3. Color: Schematic view of the structure for iron-based su- 
perconductors. Positions of pnictogen (As or P) or chalcogen (Se 
or Te) atoms are illustrated by spheres at the sites 3, 4, 7 and 
8 located either above the iron layer (scmitransparcnt sheet) de- 
picted by the plus sign inside the spheres or below the iron layer 
shown by the minus sign. Iron sites are illustrated by exam- 
ples of x 2 — y 2 (the sites 1, 6 and 9) or zx (the sites 2 and 
5) symmetries of the iron 3d orbitals, where A and B indicate 
the sublattice indices in the unit cell. The signs beside the d or- 
bitals give the rule of the local gauge employed in the present 
paper, where xy, x 2 — y 2 and 3z 2 — r 2 orbitals have the uniform 
phases, whereas zx and yz orbitals change the sign alternatingly 
depending on the A or B sublattices as is seen in the sites 2 and 
5. With this local gauge, for example, through the hybridizations 
with the pnictogen-/chalcogen-p orbitals, the transfers (includ- 
ing their signs) between the sites 1 and 2 become the same as 
those between 9 and 5 in the d model, which exemplifies the in- 
variance with the translations (±a/2,±a/2) in the ab plane, so 
that the model can have this translational symmetry apparently 
higher than the real symmetry given by the primitive translation 
(a, a) in the ab plane. 

material dependent. The x 2 — y 2 DOS has a peak at 0.3 
eV in FeTe. In FeSe, the peak is shifted to 0.7-0.8 eV 
and the pseudogap is formed. The peak position further 
shifts to 1.1 eV in LaFcAsO. The yz/zx orbitals show 
the similar trend. On the other hand, the lower peaks 
around —0.2 eV are rather pinned through the variation 
from FeTe to LaFcAsO both for x 2 — y 2 and yz / zx or- 
bitals. This leaves a pseudogap progressively in FeSe and 
then more prominently in LaFcAsO. The origin of this 
pseudogap is basically ascribed to the evolution of the h 
parameter shown in Table I. When h decreases, the hy- 
bridization between the Fe-3d and pnictogen-/chalcogen- 
p orbitals becomes appreciable and the interorbital hop- 
pings are enhanced, thus generating the splitting of the 
Fe-3cf bands into lower and upper bands. In the right 
panel of Fig. 4, the pDOS without interorbital hoppings 
is plotted. Comparison with the left panel clearly shows 
a crucial role of hybridization effect in forming the pseu- 
dogap. (we will discuss this pseudogap formation in more 
detail in Sec. 4). In the 1111 family, this large splitting 
appears to make almost a band insulator with the Fermi 
level sitting in the "gap" region between the two split 
bands for the x 2 — y 2 and yz/zx orbitals. Because of the 
incomplete band splitting, the yz/zx as well as x 2 — y 2 
orbitals contribute to the formation of the electron and 
hole pockets at the Fermi level. This semimetallic na- 
ture with small carrier density determines the relatively 
weakly correlated character of the 1111 family, together 



with the weaker effective Coulomb interaction revealed 
below. The other orbitals, 3z 2 — r 2 and xy, have clear 
gap and do not contribute to the low-energy excitations 
in all the families. In the 11 family, especially in FeTe, 
x 2 — y 2 and yz / zx bands have overall peak at the Fermi 
level even in each partial DOS and form partially filled 
bands, in the LDA picture. This large DOS at the Fermi 
level in the LDA may efficiently trigger the electron cor- 
relation effects when we consider the Coulomb interac- 
tion beyond the LDA. Indeed, it will be revealed that in 
the 11 family the effective Coulomb interaction itself is 
larger than that in the 1111 family. Combination of these 
two may lead to appreciable electron correlation effects 
when we go beyond the LDA as is known in the forma- 
tion of the Mott insulator in the half-filled band of the 
Hubbard model. We stress again that the evolution from 
the relatively weak correlation for the 1111 family to the 
strong correlation for the 11 family emerges only for the 
x 2 — y 2 and yz/zx orbitals in an orbital-selective fashion. 
While the combination of the 3z 2 — r 2 and xy orbitals is 
always "band-insulating" like and does not join in this 
physics directly, it should be noted that the formation 
of the pseudogap in x 2 — y 2 and yz/zx is caused by the 
hybridization between those orbitals and 3z 2 — r 2 and 
xy, so that it is not so trivial to derive a three-orbital 
model for x 2 — y 2 and yz/zx. 



Table II. Occupation number of the MLWFs in the d model. The 
number is normalized so that unity corresponds to half filling, 
and thus the sum of the occupancy is 6. The occupations of 
yz/zx and x 2 — y 2 are roughly close to half filling. 





xy 


yz 


3z 2 - r 2 


zx 


x 2 -y 2 


LaFcAsO 


1.14 


1.19 


1.50 


1.18 


0.99 


FeSe 


1.22 


1.06 


1.55 


1.06 


1.11 


FeTe 


1.28 


0.98 


1.53 


0.98 


1.24 



3.2 Wannier functions 

In the present work, we derive two models for each 
material. One is the d model which contains the Fe-3d 
manifold only. The other is the dp model which contains 
the P-3p/As-4p/Se-4p/Te-5p states in addition to the Fe- 
3d bands. We include 0-2p states as well in LaFePO and 
LaFeAsO, and derive the dpp model instead of the dp 
model, since they overlap with pnictogen-p band. 

Table III shows the spread (quadratic extent) of the 
MLWFs. Note that the xy axes in our convention are 
along the unit vectors of the cell containing two Fc atoms 
(Fig. 3), while they are sometimes rotated by 45 de- 
grees in other works. One striking feature is that the 
Wannier orbitals are strongly orbital dependent in the 
d model. The x 2 ~ y 2 orbital is the most extended ex- 
cept for FeTe, whereas the xy and 3z 2 — r 2 are localized. 
The anisotropy is enhanced for the 1111 family, for which 
the absolute values of the spread are large. These trends 
are understood as follows. The ten states in the d model 
contain considerable pnictogen-/chalcogen-p component. 
The hybridization between the p and Fe-d states make 
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Table III. Spread of the MLWFs (in A 2 ) defined by quadratic 
extent. The spread is the largest for the x 2 — y 2 orbital except 
for FeTc. 



LaFeAsO 



FeSe 



FeTe 





xy 


y z 


3z 2 — r 2 




x 2 — y 2 


T nPY»PO (A mnrlpl ^ 
IjcLT CIr \J \Ui IIIUQCI ) 


9 QO 


3 R7 
o.o / 


o. ry 


o.o / 


6 24 


T nTToAQO ( rl mnrlnl 


9 7^ 
^. / o 




0. 14 


3 Q1 


5 ^7 
o.o t 


R TTdt-i Afin 1 (7 Tn/~ir1ol 1 


2.68 


4.07 


2.65 


4.07 


4.22 


LiFcAs (d model) 


2.68 


3.16 


2.52 


3.16 


3.51 


FeSe (d model) 


1.82 


2.20 


1.73 


2.20 


2.60 


FeTe (d model) 


1.98 


3.51 


1.94 


3.51 


2.48 


LaFePO (dpp model) 


1.01 


1.67 


1.31 


1.67 


2.03 


LaFeAsO (dpp model) 


1.02 


1.65 


1.23 


1.65 


1.69 


BaFe2As2 (dp model) 


0.98 


1.27 


0.95 


1.27 


1.32 


LiFeAs (dp model) 


1.01 


1.17 


0.97 


1.17 


1.07 


FeSe (dp model) 


0.78 


0.86 


0.77 


0.86 


0.86 


FeTe (dp model) 


0.94 


0.97 


0.90 


0.97 


0.86 



the Wannier functions de-localized. The hybridization be- 
comes stronger in the x 2 — y 2 Wannier orbital, because 
the orbital is directed to the pnictogcn/chalcogen atoms, 
consequently the spread increases. Concerning the fam- 
ily dependence, pnictogen/chalcogen atoms approach the 
Fe plane in the order of the 11 — > 122 — > 1111 families 
(see Table I). Accordingly the 1111 family has larger hy- 
bridization effects. This trend is confirmed by comparing 
the Wannier orbitals (Fig. 5), where we can see that the 
Wannier orbital is more extended in LaFeAsO than in 
FeTe. 

Comparing the dp /dpp model with the d model, we find 
that the former has smaller spread as seen in compari- 
son of Figs. 5 and 6; In the dp/ dpp model, the Wannier 
functions are constructed from a larger number of the 
Kohn-Sham states, and therefore, the optimized Wannier 
orbitals are more localized and contain less p character 
(see Table III). They are more atomic-orbital like and 
the orbital dependence is weaker. 



LaFeAsO 



FcSc 



FeTe 




Fig. 5. Color: Isosurface of the maximally localized Wannier func- 
tion at ±0.02 bohr -3 / 2 for the Fe x 2 — y 2 orbital in the d model 
of LaFeAsO (left), FeSe (middle), and FeTe (right). This illus- 
trates how the Wannier spread shrinks from LaFeAsO to FeTe. 
The dark shaded surfaces (color in blue) indicate the positive 
isosurface at +0.02 and the light shaded surfaces (color in red) 
indicate —0.02. 




Fig. 6. Color: Isosurface of the maximally localized Wannier func- 
tion at ± 0.02 bohr -3 / 2 for the Fe-d x 2_ y 2 orbital in the dpp/dp 
model of LaFeAsO (left), FeSe (middle), and FeTe (right). The 
dark shaded surfaces (color in blue) indicate the positive isosur- 
face at +0.02 and the light shaded surfaces (color in red) indicate 
-0.02. 



3.3 Transfer integrals of d model 

Tables IV- VII show the transfer integrals, t mn (R), in 
the d model. (In the tables and this subsection, the sym- 
metry of the d orbitals is denoted as the number; 1 for 



xy, 2 for yz, 3 for 3z 



, 4 for zx, and 5 for x 



V 

orbitals.) The hopping parameters between the neighbor- 
ing Fe sites are listed in the column for (R x , R y , R z ) = 
(0.5, —0.5,0). The largest component is for (to, n)={xy, 
xy), t\\. Its value is 0.3-0.4 eV and is slightly larger in 
the 11 and 111 families than in the 1111 and 122 families, 
in accordance with shorter in-plane lattice constant. The 
next nearest hopping (R x , R y , R z ) = (1,0,0) is compa- 
rable with the nearest neighbor hopping. Especially, the 
(m, n) = (zx, zx) component, t' 4i , is larger than t\\ in the 
1111 and 122 families. The ratio t' 44 /tu, which can be 
regarded as a measure of the frustration, is larger in the 
1111 family than in the 11 family. Other transfer integrals 
associated with the yz or zx orbital are also large. These 
orbitals (in the d model) are well extended, hence hop- 
ping through the pnictogcn/chalcogen atom would con- 
tribute to the transfer integral. Longer-range hoppings 
are also nonnegligible. In fact, if we neglect transfer in- 
tegrals with less than 0.05 eV in the absolute value, the 
band width is reduced by 12 % (7 %) in LaFeAsO (FeSe). 
The transfer integrals in the c direction are sizable in 
BaFe2As2, particularly for the x 2 — y 2 orbital, indicat- 
ing three-dimensional character of the electronic states. 
They are much smaller in the 1111 family. 

Here we remark the rule to derive the transfers not 
explicitly shown in Tables IV- VII by the symmetry op- 
erations. As is shown in Fig. 3, there are two iron atoms 
(Fe-A and Fc-B) in the unit cell. Note that we can choose 
arbitrarily the phase of the transfer hoppings between Fe- 
A and Fe-B sites [i.e., t m AnB(R)] by introducing a local 
gauge to one of these irons. Let us first look at the case of 
R z =0. If we employ a common (global) coordinate for 
the irons and define a common phase of i'nAR an< ^ ^nBR 
(as was done in ref. 32), then t m AnB(R) = s x t m B n A(R) 
with s=— 1 when one and only one of to and n is yz or 
zx orbitals. Otherwise s = 1. However, as was done in 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 9 



Table IV. Transfer integrals in the d model, t mn (R x , R y , R z ), where m and n specify symmetry of d orbitals; 1 for xy, 2 for yz, 3 for 
3z 2 - r 2 , 4 for zx, and 5 for x 2 — y 2 orbitals. Symmetry operations of a y , I, and change t mn {R x , R y , R z ) to t mn (R x , —R y ,R z ), 
tmn(— Rx , — Ry , Rz) , and t mn (R y , R x , R z ) . Notice also that t mn (R) = t„ m (—R). Since the calculations are carried out using the 
4x4x4 fc-mesh, the Wannier functions are periodic with the period of (4,0,0). Because of this periodicity, the transfer integrals for 
R=[2,0,0] defined by eq.(l) are twice compared to the case in which the period of the Wannier functions is sufficiently long. We 
therefore halved the values for R=[2,0,0] in the tables, which are more appropriate as parameters for the model Hamiltonian. Units 
are given in meV. 



LaFePO 



R 



(m, n) 



[0,0,0] [i-i,o] [1,0,0] [1,-1,0] 



■i,0] [2,0,0] [0,0,f] [ 



(1, 1) 


598 


-342 


-106 


-25 


-11 


-7 


-1 


1 + 


+ 


+ 


(1,2) 





307 


172 


12 


18 








+ 




"(1,4) 


(!,3) 





-327 





-21 


-29 








-1 


+ 


+ 


(1,4) 





307 





12 


15 








-1 




"(1,2) 


(1.5) 














21 








-3 


+ 




(2,2) 


886 


231 


111 


28 


15 


12 


3 


6 + 


+ 


(4,4) 


(2,3) 





-31 





-13 


-8 








1 




"(4,3) 


(2,4) 





13 1 





-15 


-2 











+ 


(4,2) 


(2,5) 





164 





-7 


26 








3 




(4,5) 


(3,3) 


684 


130 


22 


-42 


-23 


-20 


-7 


-2 + 


+ 


+ 


(3,4) 





31 


115 


13 


10 








-1 + 




"(3,2) 


(3,5) 








-234 





■1 


11 


2 


-2 + 


+ 




(4,4) 


886 


231 


395 


28 


17 


82 


3 


1 + 


+ 


(2,2) 


(4,5) 





-164 


-77 


7 


-17 








2 + 




(2,5) 


(5,5) 


1234 


-257 


156 


-56 


33 


-27 


-2 


1 + 


+ 


+ 


LaFeAsO 




R 

(to, ro)^^~~--^ 


[0,0,0] 


[L-i-o] 


[1,0,0] 


[1,-1,0] 


[|,-|,o] 


[2,0,0] 


[o,o,f] 




1 




(1, 1) 


790 


-315 


-67 


-19 


-2 


1 


-2 


i + 


+ 


+ 


(1,2) 





253 


138 


1 


10 








o + 




"(1,4) 


(1.3) 





-301 





1 


-18 











+ 


+ 


(1,4) 





253 





1 


33 








-1 




"(1,2) 


(1,5) 














10 








-2 


+ 




(2,2) 


1099 


206 


135 


12 


9 


5 


1 


7 + 


+ 


(4,4) 


(2,3) 





-73 





-2 


-1 








2 




"(4,3) 


(2,4) 





137 





-18 


-9 








1 


+ 


(4,2) 


(2,5) 





165 





-1 


10 








3 




(4,5) 


(3,3) 


890 


72 


-13 


-38 


-15 


-18 


-6 


-2 + 


+ 


+ 


(3,4) 





73 


137 


2 


-3 








-1 + 




"(3,2) 


(3,5) 








-159 





1 


17 


3 


-3 + 


+ 




(4,4) 


1099 


206 


345 


12 


36 


70 


1 


+ 


+ 


(2,2) 


(4,5) 





-165 


19 


4 


-11 








1 + 




(2,5) 


(5,5) 


1255 


-152 


118 


-24 


30 


-28 


1 


-2 + 


+ 


+ 



ref. 18 if we attach a gauge field exp(i7r) to 4> vzB -r. and 
4>zxBR, then t mAnB (R) becomes equal to t mBnA (R) for 
all to and n. 

In addition, the transfer between two irons on the same 
sublattice t m AnA (R) = t m BnB (R) also has translational 
symmetry irrespective of the choice of this sublattice de- 
pendent gauge. Thus, if we choose the above appropriate 
gauge, the transfer integrals in these systems depend only 
on to and n and R irrespective of the sublattice of iron, 
which allows an apparent higher translational symmetry. 
The transfer constructed in this gauge gains an apparent 
square lattice symmetry without distinction of the two 
sublattice points. 

For the cases of 1111, 111, and 11 where the transla- 
tion of (0,0, 1) (I, any integer) generates the layer struc- 
ture at R z ^ identical to that at R z = shown in 
Fig. 3 in terms of the pnictogen/chalcogen positions rel- 



ative to the iron layers, whereas, in the case of 122, the 
identical layer structure is obtained only after (^,2,0) 
translation, because of the antiphase of the As positions. 
In other words, the conventional cell of 122 contains two 
iron layers because of the body center tetragonal sym- 
metry. 

Pnictogen/chalcogen positions are in-phase along the c 
axis for the 1111, 111 and 1 1 families , this higher transla- 
tional symmetry holds for general R z for the present local 
gauge exp(i7r) to 4> yzB n and 4>, xB R- Therefore, we can 
unfold the Brillouin zone (BZ). This is convenient when 
we want to solve the obtained effective lattice model nu- 
merically because the unit cell is halved with only one 
iron site contained. In Tables IV- VII, we list the transfer 
integrals for this gauge. Although the five-band model in 
this extended BZ is generally more convenient than the 
ten-band model in the original BZ for numerical model 



10 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



Table V. Transfer integrals in the d model for BaFc2As2. Notations are the same as Table IV. 



BaFe 2 As 2 



^ m ' n > — 


R 


[0,0,0] 




[1,0,0] 


[1,-1,0] 


[f.-*.o] 


[2,0,0] 


CTy I 




(i, 1) 




— 127 


—341 


—66 


— 18 


— 1 


10 


+ + 


+ 


(1, 2) 







260 


134 


—4 


-3 





+ 


"(1,4) 


(1, 3) 







-314 


3 


1 1 


-16 


-2 


+ 


+ 


(1,4) 







255 


3 


-7 


32 







"(1,2) 


(1, 5) 




-8 


1 


1 


1 


-2 





+ 




(2, 2) 




267 


209 


130 


2 


3 


6 


+ + 


(4,4) 


(2, 3) 







—89 


2 


1 


— 2 







-(4,3) 


(2,4) 







123 


2 


-27 


-18 


-1 


- + 


(4,2) 


(2.5) 







183 


-3 


-6 


11 







(4,5) 


(3,3) 




5 


45 


-16 


-38 


-18 


-15 


+ + 


+ 


(3,4) 







89 


158 


-4 


-5 





+ - 


"(3,2) 


(3,5) 










-160 





6 


17 


+ + 


— 


(4,4) 




267 


210 


356 


1 


41 


73 


+ + 


(2,2) 


(4,5) 







-178 


55 


10 


-11 





+ - 


(2,5) 


(5,5) 




363 


-133 


107 


-30 


41 


-28 


+ + 


+ 


/^>\ 
\7Yi , n ) 




[0.0,£]ab 


[0.0,£]ba 


ri i jl] 

12 ' 2 > 2a 1 


ri I 

12 ' 2 ' 2a 1 










(i, i) 




-10 


-10 


- 1 


- 1 










(1, 2) 










3 


-16 










(1, 3) 




1 


1 


2 


2 










(1,4) 










16 


-3 










(1, 5) 
























(2, 2) 




-23 


—23 


— 14 


— 14 










(2, 3) 










1 


— 7 










(2,4) 










-11 


1 










(2,5) 










-10 


-1 










(3,3) 




-86 


-86 


-37 


-37 










(3, 4 ) 










1 


-7 










(3,5) 




-111 


111 


51 


-51 










(4,4) 




-23 


-23 


-14 


-14 










(4,5) 










1 


-10 










(5,5) 




-162 


-162 


85 


85 











calculations, note that we have to fold the Brillouin zone 
again to compare the result with experiments. 

For 122, while the relation t m A n B{R) = t m B n A{R) 
holds for R z = 0, it is not satisfied for R z + 1 
because of the above mentioned antiphase of the As 
position. Namely, the BZ can not be unfolded. For 
example, i mJ 4nB(0, 0, — j^) is not necessarily equal to 
t m BnA(0, 0, — ■#-), so that we list both of them in Table 
V. 

In Tables IV-VIL we also show that how the trans- 
fer hoppings change their signs by symmetry operations. 
For example, the transfer hoppings of t mn (-^, 5, 0) is ob- 
tained by applying <j y to the transfer for R=(^, — ^,0), 
namely t mn {\, \, Q)=a y t mn {\, -\, 0), while those of 
R=(— 1,|,0) is obtained by applying / instead of a y . 
The transfer for it=(|,|,0) is obtained by apply- 
ing cry, I, and cr d to i?=(|,-i,0) as t mn (^,^,0) = 
cr y I(Tdt m n{%, — hi 0)- In the column cr^, ±(m, n) means 
that one should take it mn at the same R irrespective of 
the column index (m,n). 

3.4 Screened Coulomb interaction of d model 

Table VIII shows the onsitc screened Coulomb interac- 
tion, U mn (0), and onsite screened exchange interaction, 



J„(0), in the d model. As has been reported previously 
for the 1111 family, 32 ) U depends strongly on the orbital 
by the amount as large as >1 eV. The orbital dependence 
is seen in other families as well. The values of U decrease 
and the anisotropy is enhanced as the Wannier orbital is 
extended. 

Table IX shows the average of the diagonal (m = n) 
terms of U. It is the smallest in the 1111 family. The 
value increases slightly by ^0.3 eV in the BaFe2As2, and 
is substantially larger in the 111 and 11 families. The 
relative interaction strength U ft in the d model is nearly 
8 for LaFcPO, 9 for LaFcAsO, and 14 for FeSe. There- 
fore FeSe is substantially more strongly correlated. Here, 
U is defined as the average of the diagonal intraorbital 
onsite interactions. When we specify the largest transfer 
integral i ma xri among pairs between an orbital n at an 
iron site and any 3d orbitals at its nearest-neighbor iron 
site, t is defined as the average of t maxn over the five or- 
bitals; i = J2n=i ^maxn/5- For example, t for LaFeAsO is 
defined by (0.32+0.25+0.30+0.25+0. 17)/5=0.26. 

Small (large) U in the 1111 (11) family originates from 
the following two factors: One is the MLWF basis for 
which U is defined. The MLWF is more extended (local- 
ized) in the 1111 (11) family, thereby the bare interaction 
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Tabic VI. Transfer integrals in the d model for LiFcAs. Notations are the same as Table IV. 



LiFeAs 



(m, 



R 


[0,0,01 


[5.-5.0] 

L 2 ' 2 ' J 


[1,0,0] 


[1,-1,0] 

L ' 'J 


[1.-5.0] 

L 2 ' 2 ' J 


[2,0,0] 

L ' ' J 


[0,0, £] 


L2 ' 2 ' aJ y 


I 




(i, i) 


-180 


-409 


-26 


-57 


1 


25 


-32 


9 + 


+ 


+ 


(1,2) 





290 


144 


-3 


-15 








-10 + 


- 


"(1,4) 


(1,3) 





-346 





35 


-8 








17 


+ 


+ 


(1,4) 





290 





-3 


56 








-3 


- 


"(1,2) 


(1,5) 














-19 








-7 


+ 


- 


(2,2) 


270 


229 


142 


-32 


-21 


-9 


11 


39 + 


+ 


(4,4) 


(2,3) 





-127 





11 


4 








6 




"(4,3) 


(2,4) 





167 





-42 


-30 








15 


+ 


(4,2) 


(2.5) 





198 





-9 


-3 








18 




(4,5) 


(3,3) 


-133 


-26 


-52 


-13 


5 


-19 


-65 


-31 + 


+ 


+ 


(3,4) 





127 


188 


-11 


-18 








-35 + 




"(3,2) 


(3,5) 








-79 





-5 


17 


46 


-42 + 


+ 




(4,4) 


271 


230 


ill 


-32 


52 


93 


11 


5 + 


+ 


(2,2) 


(4.5) 





-198 


124 


9 


-19 








30 + 




(2,5) 


(5,5) 


137 


-7 


81 


1 


35 


-37 


58 


-51 + 


+ 


+ 



v becomes small (large) . The other factor is the strength 
of the screening. As is shown in Table IX, the ratio of the 
screened Coulomb interaction U to the bare interaction v 
is smaller in the 1111 family than in the 11 family. This 
can be understood as follows. Firstly, there is a larger 
number of bands near the Fe-d band, such as occupied 
O-p band and dense unoccupied states above the Fe-d 
band. They contribute to screening in the 1111 systems. 
Secondly, the energy levels of the pnictogcn-p states in 
the 1111 systems are shallower than the Se-p level in 
FeSe. It also enhances the screening. Hybridization be- 
tween the d and pnictogen/chalcogen orbitals would also 
affect the screening effect. As the hybridization becomes 
stronger, the transition matrix element would become 
larger, which makes screening more effective. 

There is little difference in U between LaFcPO and 
LaFeAsO, although U is slightly larger in LaFcAsO. On 
the other hand, in the 11 system, FeSe is significantly 
more correlated than FeTe. The values of U, v, and U/v 
are 24 %, 4 %, and 20 % larger in FeSe than in FcTc, 
respectively. This means that the screening effect in W r , 
not the size of the MLWFs, is curial for the difference. 
The chalcogen-p states are shallower in FeTe than in 
FeSe. Transition energies between the p states and unoc- 
cupied states are thus smaller, which results in stronger 
screening in FeTe. Also, the p bands are entangled with 
the Fe-d bands. This would change the transition matrix 
elements and leads to the stronger screening. 

The p-d hybridization delocalizes the Fe-3d yz and zx 
Wannier orbitals in FeTe. In contrast to the other com- 
pounds, their spreads are the largest among the five or- 
bitals. It is larger by 60 % than that in FeSe, and close to 
that of LaFePO (see Table III). U is accordingly smaller 
than other orbitals. 

We note that the value of U in LaFeAsO is somewhat 
smaller than that in the previous works. 61, 62 -* We found 
that energy levels of the unoccupied states are sensitive 
to the number of basis functions taken into account in 
the FP-LMTO band calculations. Inclusion of more ba- 



sis functions decreases U slightly compared to the pub- 
lished data. On top of that, we found that the presence 
of the localized La-4/ states in the 1111 family quantita- 
tively affects the resulting U . In the quantitative aspect, 
in LDA, the La-4/ levels are expected to be too low, 
thus leading to artificially smaller U . In order to check 
this effect, wc performed the cRPA calculation starting 
from the LDA+J7 solution where the input-!/ of 1 Ry is 
imposed on La-4/ states (this input- U pushed up La-4/ 
level by ~ 5 eV). The resulting effective onsite interaction 
U for the Fe-3c? orbitals was found to increase roughly by 
0.2 eV. We will come back to this point and examine it 
in more details later. 

The average of the exchange energies is 0.4 eV in the 
1111, 122, 111 families. In the 11 family, the value is 
slightly larger and is 0.5 eV. This trend is the same as 
the Coulomb interaction, though the family dependence 
is much weaker. As shown previously for the 3d transition 
metals, 66 - 1 screening effects arc small for the exchange en- 
ergy. The value is reduced only by ^0.1 eV compared to 
the bare exchange interaction in all the compounds, in 
sharp contrast with the diagonal Coulomb interaction, 
where the reduction due to screening is a factor of 4-6. 
The family dependence primarily comes from the differ- 
ence in the spatial extent of the MLWFs rather than the 
strength of screening. 

Interaction between the nearest neighbor Fe sites, V, 
is shown in Table X. The value of V is 0.6-1.1 eV, and the 
orbital dependence of V is weak in contrast to the onsite 
U. The crucial factor for V is not the spatial extent of 
the MLWFs but the distance between the sites. Longer- 
range interactions are found to decay as 1/ar, where r is 
the distance between the Wannier centers and a=0.5-0.6 
(eV-A)- 1 in the 1111 and 122 families, and a=0.3.-0.4 
(eV-A) -1 in the 11 and 111 families. 

The presented results for the d model are summarized 
in Fig. 7. From this plot, the systematic change of the 
screened interaction and the geometrical frustration mea- 
sured from the amplitude of the next-nearest-neighbor to 
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Table VII. Transfer integrals in the d model for FeSe and FcTe. Notations are the same as Table IV. 



FcSc 



R 



(m, n) 



[0,0,0] [|-|,0] [1,0,0] [1,-1,0] 



■1,0] [2,0,0] [0,0,f] [ 



(1, 1) 


854 


-410 


-69 


-10 


4 


11 


-24 


7 + 


+ 


+ 


(1,2) 





272 


131 


-9 


-6 








-8 + 


- 


"(1,4) 


(1,3) 





-347 





22 


-7 








11 


+ 


+ 


(1,4) 





272 





-9 


23 








-2 


- 


"(1,2) 


(1,5) 














-8 








-4 


+ 




(2,2) 


1418 


199 


126 


-16 


-8 


-6 


9 


30 + 


+ 


(4,4) 


(2,3) 





-120 





7 


4 








10 




"(4,3) 


(2,4) 





127 





-24 


-18 








11 


+ 


(4,2) 


(2.5) 





223 








-4 








20 




(4,5) 


(3,3) 


980 


-3 


-18 


-16 


-3 


-15 


-22 


-8 + 


+ 


+ 


(3,4) 





120 


196 


-7 


-12 








-10 + 




"(3,2) 


(3,5) 








-115 





-3 


10 


-7 


-6 + 


+ 




(4,4) 


1418 


190 


348 


-16 


15 


59 


9 


1 + 


+ 


(2,2) 


(4.5) 





-223 


82 





-13 








7 + 




(2,5) 


(5,5) 


1335 


-57 


92 


-1 


19 


-24 


-29 


5 + 


+ 


+ 


FeTe 




(m, n) 


[0,0,0] 


[i-Lo] 


[1,0,0] 


[1,-1,0] 


[§,-*.o] 


[2,0,0] 


[o,o,f] 


[ 2 ' 2 ' a ] a V 


I 


0",; 


(1,1) 


163 


-392 


-2 


-14 


2 


8 


-39 


10 + 


+ 


+ 


(1,2) 





228 


06 


-9 


-13 








-15 + 




"(1,4) 


(1,3) 





-341 





33 


-6 








21 


+ 


+ 


(1,4) 





228 





-9 


39 








1 




"(1,2) 


(1, 5) 














-15 











+ 




(2,2) 


774 


164 


99 


-33 


-17 


-15 


10 


37 + 


+ 


(4,4) 


(2,3) 





-130 





16 


8 








15 




"(4,3) 


(2,4) 





107 





-31 


-34 








17 


+ 


(4,2) 


(2,5) 





181 





-5 


-9 








18 




(4,5) 


(3,3) 


189 


-87 


-68 


4 


19 


-18 


-47 


-12 + 


+ 


+ 


(3,4) 





130 


200 


-16 


-11 








-30 + 




"(3,2) 


(3,5) 








-23 





-13 


13 


15 


-20 + 


+ 




(4,4) 


774 


164 


348 


-32 


17 


80 


10 


11 + 


+ 


(2,2) 


(4,5) 





-184 


111 


5 


-26 








25 + 




(2,5) 


(5,5) 


466 


85 


40 


9 


10 


-16 


17 


-25 + 


+ 


+ 



the nearest-neighbor transfers are more visible with the 
evolution from LaFcPO, LaFeAsO, BaFc2As2, LiFeAs to 
FeSe and FeTe. We again emphasize that the family de- 
pendence of U/i comes from U not from i. 

3.5 Transfer integrals of dp/dpp model 

The effective parameters are very different in the dpp 
(in LaFePO and LaFeAsO) or dp (in other materials) 
model compared to those in the d model. The change 
in the spread and shape of the MLWFs (compare Figs. 5 
and 6) results in the change in the transfer integrals. The 
largest transfer integral between the nearest-neighbor d 
orbitals, t max , is the d xy -d xy hopping. Its value is compa- 
rable or even larger in the dpp /dp model; 0.316 eV, 0.403 
cV, 0.367 eV for LaFeAsO, FeSe and FeTe, respectively, 
against 0.315 eV, 0.410 eV, 0.392 eV in the d model. 
Other nearest d-d transfer integrals alter significantly, 
because the change in the MLWFs is significant. Most 
of them are smaller in magnitude in the dpp /dp model. 
However, the d x 2_ y 2-d x 2_ y 2 transfer integral gets larger. 
The value is 0.243 eV, 0.272 cV, 0.308 cV for LaFeAsO, 
FeSe, FeTe, respectively in the dpp/ dp model, while it is 



0.152 eV, 0.057 eV, 0.085 eV in the d model. The next- 
nearest-neighbor d-d transfer integrals are small and even 
the largest value is ~0.1 eV in all the materials, to be 
compared with 0.3-0.4 eV in the d model, clearly indi- 
cating that the next-neighbor transfers in the d model is 
mediated by the pnictogen-/chalcogen-p orbitals. In fact, 
in the dp/dpp model, the transfer integrals between the 
Fe-d and pnictogen- / chalcogen-p orbitals are larger than 
the d-d transfer integrals; the largest one is the d zx -p x 
transfer integral, the value of which is 0.734 cV, 0.895 
eV, and 0.711 cV for LaFeAsO, FeSe and FeTe, respec- 
tively. It should be noted here that, in comparison of the 
three compounds, the value is the largest in FcSc. This 
is unexpected; the spread of the d zx orbital is the largest 
and the distance between Fe and pnictogen/chalcogen 
atoms is the shortest in LaFeAsO. Therefore, naively we 
expect that the the transfer integral would be the largest 
in LaFeAsO, which turned out not to be the case. Con- 
cerning the p-p transfer integrals between the nearest 
pnictogen/chalcogen sites, the largest value is for the p z - 
p z transfer, which is 0.306 eV, 0.301 eV, 0.319 eV in 
LaFeAsO, FeSe, and FeTe, respectively. The p x -p x (p y - 
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Table VIII. Effective on-site Coulomb ( U) /exchange (J) interactions between two electrons on the same iron site in the d model for all 
the combinations of iron-3d orbitals (in eV). 



LaFcPO 






U 












J 








xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 




xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 


X if 


2.98 


1.80 


1.78 


1.80 


1.77 


XL) 




0.45 


0.54 


0.45 


0.20 


u" 


1.80 


2.42 


1.97 


1.64 


1.46 




0.45 




0.32 


0.36 


0.31 


3^2 _ r 2 


1.78 


1.97 


2.81 


1.97 


1.46 


3^2 _ r 2 


0.54 


0.32 




0.32 


0.37 


zx 


1.80 


1.64 


1.97 


2.42 


1.46 


zx 


0.45 


0.36 


0.32 




0.31 


x 2 — y 2 


1.77 


1.46 


1.46 


1.46 


1.68 


x 2 — y 2 


0.20 


0.31 


0.37 


0.31 




LaFcAsO 






U 












J 








xy 


yz 


q ~2 '2 
OZ — T 


Z X 


2 2 
x — y 




x y 


yz 


OZ — T 


ZX 


^2 ,,2 

x — y 


xy 


3.03 


1.80 


1.78 


1.80 


1.91 


xy 




0.46 


0.57 


0.46 


0.23 


yz 


1.80 


2.43 


1.97 


1.62 


1.52 


yz 


0.46 




0.33 


0.37 


0.35 


3z 2 - r 2 


1.78 


1.97 


2.84 


1.97 


1.51 


3z 2 - r 2 


0.57 


0.33 




0.33 


0.42 


zx 


1.80 


1.62 


1.97 


2.43 


1.52 


Z X 


0.46 


0.37 


0.33 




0.35 


x 2 — y 2 


1.91 


1.52 


1.51 


1.52 


1.91 


x 2 — y 2 


0.23 


0.35 


0.42 


0.35 




BaFeAs 






IJ 












J 


2 






xy 


yz 


^ 9 9 

3z 2 - r 2 


zx 


9 9 

x 2 - y 2 




xy 




^.9 9 

3z 2 — r 2. 


zx 


9 9 

x — y 


xy 


3.18 


1.94 


1.99 


1.94 


2.16 


xy 




0.48 


0.60 


0.48 


0.26 


yz 


1.94 


2.64 


2.21 


1.77 


1.72 


yz 


0.48 




0.36 


0.40 


0.41 


3z 2 - r 2 


1.99 


2.21 


3.28 


2.21 


1.77 


3z 2 - r 2 


0.60 


0.36 




0.36 


0.50 


zx 


1 O/l 

i.y4 


1 77 




O P.A 


1 70 


Z X 


U.4o 


n a n 
U.4U 


U.oD 




n a 1 
U.41 


x 2 — y 2 


2.16 


1.72 


1.77 


1.72 


2.29 


x 2 — y 2 


0.26 


0.41 


0.50 


0.41 




LiFeAs 






u 












J 








xy 


yz 


o 9 9 

3z 2 — r 2 


zx 


9 9 

x 2 — y 2 




xy 


yz 


n 9 9 

3z z — r 2. 


zx 


9 9 

x 2 - y 2 


xy 


3.39 


2.23 


2.27 


2.23 


2.54 


xy 




0.47 


0.60 


0.47 


0.28 




2.23 


2.96 


2.52 


2.08 


2.11 


yz 


0.47 




0.36 


0.39 


0.44 


3z 2 - r 2 


2.27 


2.52 


3.58 


2.52 


2.15 


3z 2 - r 2 


0.60 


0.36 




0.36 


0.54 


zx 


2.23 


2.08 


2.52 


2.96 


2.11 


zx 


0.47 


0.39 


0.36 




0.44 


X 2 — V 2 

■' y 


2.54 


2.11 


2.15 


2.11 


2.85 


X 2 — V 2 


0.28 


0.44 


0.54 


0.44 




FcSc 






u 












J 








xy 


yz 


r. 9 9 

3z 2 - r 1 


zx 


9 9 

x z — y z 




xy 


yz 


f-. 9 9 

3z 2 — r 2 


zx 


9 9 

x — y 


xy 


4.51 


3.19 


3.20 


3.19 


3.49 


xy 




0.57 


0.69 


0.57 


0.32 


yz 


3.19 


4.11 


3.52 


3.02 


2.98 


yz 


0.57 




0.42 


0.48 


0.53 


3z 2 - r 2 


3.20 


3.52 


4.67 


3.52 


3.00 


3z 2 - r 2 


0.69 


0.42 




0.42 


0.62 




3.19 


3.02 


3.52 


4.11 


2.98 




0.57 


0.48 


0.42 




0.53 


x 2 - y 2 


3.49 


2.98 


3.00 


2.98 


3.78 


x 2 - y 2 


0.32 


0.53 


0.62 


0.53 




FeTe 






U 












J 








xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 




xy 


yz 


3z 2 - r 2 


zx 


x 2 -y 2 


xy 


3.84 


2.34 


2.50 


2.34 


3.04 


xy 




0.49 


0.68 


0.49 


0.34 


yz 


2.34 


2.88 


2.56 


2.03 


2.29 


yz 


0.49 




0.37 


0.37 


0.49 


3z 2 - r 2 


2.50 


2.56 


3.84 


2.57 


2.44 


3z 2 - r 2 


0.68 


0.37 




0.37 


0.66 


zx 


2.34 


2.03 


2.57 


2.88 


2.29 


zx 


0.49 


0.37 


0.37 




0.49 


x 2 — y 2 


3.04 


2.29 


2.44 


2.29 


3.59 


x 2 — y 2 


0.34 


0.49 


0.66 


0.49 





Table IX. Average of the diagonal terms of the screened Coulomb interaction, U, and that of the bare Coulomb interaction v. As a 
measure for the degree of screening, the value of U /v is also shown. 





d model 


dp/dpp model 




U (eV) 


v (eV) 


U/v 


U (eV) 


v (eV) 


U/v 


LaFcPO 


2.47 


14.15 


0.174 


4.13 


18.96 


0.218 


LaFeAsO 


2.53 


14.85 


0.171 


4.23 


19.46 


0.217 


BaFe 2 As 2 


2.80 


15.59 


0.180 


5.24 


20.38 


0.257 


LiFeAs 


3.15 


15.82 


0.199 


5.94 


20.35 


0.292 


FeSe 


4.24 


17.53 


0.242 


7.21 


21.37 


0.337 


FeTe 


3.41 


16.89 


0.202 


6.25 


20.90 


0.299 
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Fig. 7. Material dependence of parameters in d model. The average of the onsite effective Coulomb interactions (U), the average of the 
offsite effective Coulomb interactions between the neighboring Fe sites (V), the average of the onsite effective exchange interactions (J), 
the maximum value of the transfer integrals between the neighboring Fe sites [tn = t n (1/2, — 1/2, 0)] and between the next-nearest 
neighbor [t' i4 = t44 (1,0,0)], U/i, and t' 4i /tn arc compared. The subscripts of tn and t' i4 are orbital indices; 1 for xy and 4 for zx. i 
is the orbital average of the largest nearest d-d transfer integrals. 



p y ) hopping is 0.255 eV, 0.203 eV and 0.164 eV, and the 
Px-Pz (Py-Pz) hopping is 0.279 eV, 0.272 eV and 0.225 
eV in each material. 

3.6 Screened Coulomb interaction ofdp/dpp model 

Table XI presents the U and J matrices in the dp/dpp 
model. The value of U is large compared to the d model 
for the following two reasons. 62, 66 ^ Firstly, there are more 
states kept in the model, so that more screening processes 
are taken away from W r . Secondly, the MLWFs are more 
localized, since they are optimized using more states. In 



addition, the MLWFs in the dpp/dp model is less orbital 
dependent. Consequently, the orbital dependence in U is 
much weaker than the d model. The exchange value is 
somewhat larger than that in the d model by ~0.2 eV. 

Examining the family dependence, we find that U is 
substantially larger in the 11 family than that in the 
1111 family: The value of U/i (U/W) is 20 (0.48), 20 
(0.58), 30 (0.91), 26 (0.77) in LaFcPO, LaFcAsO, FeSe, 
FcTc, respectively Here, U is the average of the diago- 
nal onsite interactions between the d orbitals, and W is 
the width of the d band. We take the definition of i as 
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Table X. Effective Coulomb interaction between the neighboring sites, V, in the d model (in eV). 



LaFePO 


xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 


LaFeAsO 


xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 


xy 


0.77 


0.74 


0.74 


0.74 


0.74 


xy 


0.73 


0.70 


0.69 


0.70 


0.71 


vz 


0.74 


0.72 


0.72 


0.73 


0.71 


u z 


0.70 


0.68 


0.68 


0.69 


0.69 


3z 2 - r 2 


0.74 


0.72 


0.71 


0.72 


0.71 


3z 2 - r 2 


0.69 


0.68 


0.67 


0.68 


0.68 


zx 


0.74 


0.73 


0.72 


0.72 


0.71 


zx 


0.70 


0.69 


0.68 


0.68 


0.69 


x 2 — y 2 


0.74 


0.71 


0.71 


0.71 


0.71 


x 2 - y 2 


0.71 


0.69 


0.68 


0.69 


0.71 


BaFe2As2 


xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 


LiFcAs 


xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 


xy 


0.66 


0.63 


0.64 


0.63 


0.66 


xy 


0.94 


0.92 


0.91 


0.92 


0.93 


yz 


0.63 


0.61 


0.62 


0.63 


0.63 


yz 


0.92 


0.91 


0.90 


0.92 


0.92 


3z 2 - r 2 


0.64 


0.62 


0.62 


0.62 


0.63 


3z 2 - r 2 


0.91 


0.90 


0.89 


0.90 


0.91 


z X 


0.63 


0.63 


0.62 


0.61 


0.63 


zx 


0.92 


0.92 


0.90 


0.91 


0.92 


x 2 — y 2 


0.66 


0.63 


0.63 


0.63 


0.66 


x 2 — y 2 


0.93 


0.92 


0.91 


0.92 


0.95 


FeSe 


xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 


FeTe 


xy 


yz 


3z 2 - r 2 


zx 


x 2 — y 2 


xy 


1.11 


1.09 


1.08 


1.09 


1.11 


xy 


0.94 


0.92 


0.92 


0.92 


0.94 


yz 


1.09 


1.08 


1.07 


1.10 


1.09 


yz 


0.92 


0.91 


0.91 


0.92 


0.91 


3z 2 - r 2 


1.08 


1.07 


1.06 


1.07 


1.08 


3z 2 - r 2 


0.92 


0.91 


0.91 


0.91 


0.91 


zx 


1.09 


1.10 


1.07 


1.08 


1.09 


zx 


0.92 


0.92 


0.91 


0.91 


0.91 


x 2 — y 2 


1.11 


1.09 


1.08 


1.09 


1.11 


x 2 — y 2 


0.94 


0.91 


0.91 


0.91 


0.94 



the same as the case of the d model; t is defined by the 
largest transfers between the nearest-neighbor iron sites 
averaged over the 3d orbitals. As is shown in Table IX, 
the bare Coulomb interaction v differs only by 10 % in 
the two families, while U/v of the 1111 family is smaller 
by a factor more than 25 % than that of the 1111 fam- 
ily. Therefore, the family dependence of U prominent in 
Table IX is determined primarily by the difference in the 
screening while the spatial extent of MLWFs plays an 
only supplementary role. 

The effective onsite Coulomb interaction of the 
pnictogen-/chalcogen-p orbitals, U p , are also large (not 
shown in the table). It is the smallest in LaFePO, still the 
average of the diagonal terms is 2.5 eV. The value is the 
largest for FeSe at 4.7 eV. Interaction between the Fe-ci 
and neighboring-pnictogen-p orbitals, U pd , is not negli- 
gible as well, ranging from 1.2 eV (LaFeAsO) to 1.7 eV 
(FeSe). 

3. 1 Comparison between FP-LMTO and pseudopoten- 
tial calculations 
For the critical check of the reliability and convergence 
of the derived parameters, we compared the FP-LMTO 
results with results obtained from ab initio pseudopoten- 
tial calculations with the plane-wave basis. Comparisons 
were made for onsite U and J parameters of LaFeAsO 
and FeSe. In LaFeAsO, the pseudopotcntial results (top 
two left 5x5 matrices in Table XII) give somewhat larger 
values by 0.1-0.2 eV than the FP-LMTO values listed in 
Table VIII. This is because the present pseudopotential 
calculation does not include La-/ states and therefore 
the screenings from the La-/ states were completely ne- 
glected in the RPA calculations. On the other hand, the 
FP-LMTO results include the La- /-screening effects and 
thus give reasonably smaller values than the pseudopo- 
tential ones. It should be noted here that, in general, 
LDA tends to underestimate the energy difference be- 
tween the level of the localized / state and the Fermi 
level. In the present calculation, this level is located 
around 3 eV above the Fermi level [see Fig. 1 (a)], which 



may be too low. In order to analyze the effect of the 
La-/ level on the derived parameters, we performed con- 
strained RPA calculations in LDA+J7 formalism with the 
FP-LMTO implementation, where we employed 1 Ry for 
input-t/ on the La-/ orbitals, which pushed the / level 
up by nearly 5 eV from the original position. The re- 
sultant U and J parameters are shown in the top two 
right 5x5 matrices in Table XII, from which we see that, 
with this input-[7, LDA+J7 gives values very similar to 
the pseudopotential results and consistently suggests this 
possible small correction (~ 0.2 eV). 

In contrast, in FeSe, there exists no ambiguity due to 
the La-/ state; the pseudopotential results should agree 
with the original FP-LMTO results without resorting to 
LDA+U. Comparisons are shown in the bottom two part 
in Table XII. We see an excellent agreement between the 
two results, which confirms that the constrained RPA re- 
sults based on the MLWFs are neither affected by details 
of treatments of core electrons, nor by the basis choices 
(i.e. either plane waves or FP-LMTO), if the cutoff radius 
in the pseudopotential is taken small enough, and reason- 
ably large number of unoccupied states are included in 
the cRPA calculation. 

We here note on the convergence of the cRPA. In the 
cRPA, the number of bands participating in the screening 
should be taken sufficiently large in the part away from 
the Fermi level, so that the polarization calculation con- 
verges. In the present pseudopotential calculations, the 
total number of conduction bands, Nb a „d, is taken up to 
130 together with the number of the valence bands, 34 
for LaFeAsO and 18 for FeSe. With this choice of the 
numbers, we infer that the underestimate error of the 
polarization may lead to an overestimate of the diagonal 
part of the screened interaction with the amount at most 
0.1 eV. (For example, for U xy of LaFeAsO, the value is 
3.43 eV for N band = 50, 3.21 eV for N band = 70, and 3.14 
eV for Nband — 130.) If one wishes to reach better accu- 
racy, one needs to take larger number of bands, while for 
the present purpose with the accuracy of the order of 0.1 
eV, the present choice may be sufficient. 
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Tabic XI. Effective Coulomb (U) / exchange (J) interaction between two electrons at 3d orbitals on the same iron site in the dpp (for 
LaFePO and LaFcAsO) or dp (in other compounds) model (in eV). 
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Finally, we briefly comment an effect of replacement 
of LDA by GGA on the derived parameters. We found 
that such a replacement hardly affects the resultant pa- 
rameters; GGA gives slightly larger values than those 
obtained with LDA. The difference is less than 0.05 eV 
at the maximum and thus the choice of LDA or GGA 
would not be essential in the derivation of the model 
parameters. 

4. Summary and Discussion 

In this paper, we have derived effective low-energy 
models of iron-based superconductors, by applying the 
three-stage scheme including the downfolding to elim- 
inate the higher-energy excitation channels and to re- 
tain the low-energy (target bands) degrees of freedom 
near the Fermi level. The models derived in this paper 
are those for LaFePO, and LaFeAsO belonging to the 



1111 family, BaFe 2 As 2 (122) and LiFeAs (111) as well 
as FeSe and FeTe belonging to the 11 family. For each 
compound, two different types of low-energy models are 
presented: One, called the d model, is constructed from 
the bands which contain mainly Fe-3c? bands only. In this 
case, the model contains 10 orbitals per unit cell con- 
taining two iron sites. In the other model called the dp 
model, As-4p/P-3p/Se-4p/Te-5p orbitals arc additionally 
included in the low-energy models. In case of LaFePO 
and LaFeAsO, 0-2p orbitals are included as well, which 
constitute the dpp models. The downfolding procedure 
starts from the LDA calculation results of the global 
band structures. The MLWFs are constructed for the tar- 
get low-energy bands. The screened Coulomb interaction 
between two electrons in the target bands are derived by 
the constrained RPA, where the screening effects arising 
from the eliminated bands are included. The effective 
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Tabic XII. Effective onsite Coulomb/exchange interactions in the d model obtained with ab initio pseudopotential (PP) calculations, 
compared with the results by the FP-LMTO. The FP-LMTO calculation for LaFcAsO is done using the LDA+C/ band structure with 
input-t/ of 1 Ry for the La-4/ orbitals, leading to a shift of the La-4/ level by ~ 5 eV from the original position. 
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low-energy models are derived from the two independent 
LDA calculations, one based on the FP-LMTO and the 
other based on the pseudopotential with the plane wave 
basis. These two methods give very good agreements 
and the resulting model parameters do not depend on 
the choice of the basis function in the LDA calculations, 
which assures the reliability of the downfolding scheme 
presented here. We also checked that it neither depends 
on the choice of LDA or GGA. 

In addition, in the constrained calculations, we have 
applied a recently developed technique to the disentan- 
glement of bands 67 ) in the procedure of the downfolding 
to disconnect the target band from other entangled band 
structure. This is particularly helpful for FeTe, where Fc- 
3d bands are entangled with the Te-p bands. 

The band structures of the six compounds share global 
similarity. All of them at the Fermi level show small elec- 
tron pockets around the M point and small hole pockets 
around the T point. However, there exist nonnegligible 
differences and dependence of the derived model param- 
eters on the compounds and families. We start from the 
overall family dependence that does not depend on the 
choice of the d or dp/dpp models. One important ori- 
gin of the family dependence of the model parameters 
is the variation of the distance h between the pnictogen 
or chalcogen atoms and the Fe layer, as listed in Table 
I. In the ascending order, the distance becomes longer 
from the 1111, 111, 122 to 11 families. The distance is 
the shortest for the LaFePO and sequentially increases in 
the direction from LaFeAsO, LiFeAs, BaFe2As2 to FeSe 



and FeTe. Then, in this order, the compounds progres- 
sively lose covalent character of the chemical bonding be- 
tween Fe and the chalcogcn/pnictogen elements and gain 
stronger ionic character. In this view, FeTe is expected 
to have the largest ionic character. On top of that, elec- 
tronegativity is also important for analyses of the chem- 
ical bonding. The value of electronegativity for each el- 
ement is 1.8 for Fe, 2.1 for P, 2.0 for As, 2.4 for Se, and 
2.1 for Te. So, in this view, the bonding of FeSe has the 
largest ionic character. By contrast to FeSe and FeTe, the 
1111 family has relatively strong hybridization between 
Fe-3d and pnictogen-p orbitals. Indeed, the difference of 
the distance h, as well as the electronegativity, explains 
many of the family dependence as revealed below in a 
unified fashion. 

The 11 family, in particular, FeTe has more entangled 
band structure of the Fe-3d and chalcogen-p bands, lead- 
ing to the smearing of the pseudogap structure within 
the 3d bands above the Fermi level in contrast to the 
prominent pseudogap (located at ~ 0.5 eV higher than 
the Fermi level) in the 1111 family. In terms of the d 
model, the origin of the pseudogap structure can be un- 
derstood as follows. Let us first see why the 3z 2 — r 2 and 
xy bands have a gap around the Fermi level commonly 
for the four families (see Fig. 4). The energy scale (band 
width) of the xy band is the largest among the 3d bands, 
since the xy orbital points toward neighboring Fe sites. 
As is well appreciated, any two-dimensional single-orbital 
tight-binding models do not make gaps in the DOS unless 
the translational symmetry is broken. On the other hand, 
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as mentioned in the previous section, all the families of 
the iron compounds have a high translational symmetry 
(for the appropriate choice of the sublattice dependent 
gauge). Therefore, the gap in the xy band is not the 
single-band origin of the sole xy orbital. Instead, the ori- 
gin is explicable only with the hybridization between xy 
and other 3d orbitals. As we can see in Tables IV- VII, 
the nearest-neighbor transfer integrals [t mn (^, |,0)] be- 
tween m = xy and n = yz, zx, 3z 2 — r 2 arc appreciable; 
they are all ~ 300 mcV. It should be noted here that, 
even the 3z 2 — r 2 orbital perpendicular to the iron layer 
has rather large hybridization; t xy3z 2_ r 2 (i, ^,0) ~ 300 
meV. Through this hybridization, the xy and 3z 2 — r 2 
bands tend to make a clear gap around the Fermi level. 

For the x 2 — y 2 orbital, the nearest t xyx 2_ y 2 (i i, 0) is 
always small, but the next-nearest t x 2_ y 2 - iz 2_ r 2 (1, 0, 0) 
can be large for systems with smaller h. For example, 
i x 2_j / 2 i3z 2_ r 2 (1,0,0) is just -23 meV for FcTc but -234 
meV for LaFePO. This is because the Wannier spread of 
the x 2 — y 2 orbital becomes large for smaller h (see Tabic 
III) . The pscudogap structure in the x 2 — y 2 band thus 
appears for small-/i system. We note that the two-peak 
structure is observed in the x 2 — y 2 DOS even without 
interorbital hybridization, as shown in the right panel 
of Fig. 4. However, the energy level of the lower peak is 
shifted up in that case, thus the pseudogap is not formed 
at the Fermi level. Also, the energy separation between 
the two peaks is too small. Analysis using the dpp/dp 
model even more clearly reveals that the psudogap in- 
deed originates from the hybridization gap, indicating 
the significance of the interorbital hybridization : With- 
out the hybridization between the x 2 — y 2 and the anion 
p orbitals, which is included in the x 2 — y 2 Wannier in 
the d model, the upper peak and dip structure itself dis- 
appears. 

A similar argument can be applied to yz/zx and xy: 
When h is small, the hybridization between yz/zx and xy 
becomes strong and the yz/zx bands have a gap around 
the Fermi level. In the discussion for the yz/zx band, hy- 
bridizations of the nearest t yz / zx3z 2_ r 2 and ty Z /zx,x 2 -y 2 
are also important, especially for the 11 family. In prin- 
ciples, these transfers must vanish from the mirror-plane 
symmetry of these orbitals, but, in fact, because of the 
symmetry braking due to the presence of the anion-p 
component in the Wannier functions, the transfers ex- 
hibit finite values. Now, when h is large as in the case of 
FeTe, the symmetry tends to be largely broken so that 
tyz/zx,3z 2 -r 2 an d t yz /zx,x 2 -y 2 become large. These trans- 
fers can induce a hybridization gap in the yz/zx bands, 
but a more important point is that those cause a new 
gap formation at much lower energy than the Fermi level 
[around —1 eV, see Fig. 4(c)]. As a result, this lower gap 
formation leads to closure of the pseudogap around the 
Fermi level. 

Consequently, the 1111 family with smaller h is closer 
to the band insulator due to the band splitting formed 
by the hybridization between the 3d orbitals, regarded 
as a system with small carrier number at semimetallic 
electron and hole pockets. On the other hand, the 11 
family has a large density of states over the Fermi 
level with fewer band splitting, which may make the 



electron correlation effect more efficient, because of 
larger effective carrier number, as in the formation of 
the Mott insulator in the middle of the band in the 
Hubbard model. We note that this remarkable trend 
from the semimetal to the "half-filled" band emerges 
particularly for yz/zx and x 2 — y 2 orbitals in an orbital 
selective way. An orbital selective crossover from the 
band-insulating to Mott physics may specifically emerge 
in these yz/zx and x 2 — y 2 orbitals, while the other two 
orbitals keep more or less the band insulating character 
through this family variation. 

The dispersion of the Fe-3<i band of the 11 family has 
larger three dimensionality. This larger dispersion in the 
11 family is mainly due to the lack of the block layer 
including La and O, which makes the three-dimensional 
overlap of Fe-3<i orbital through the chalcogcn-p orbital 
larger. It appears to be triggered by the larger h of 
the 11 family as well. The larger h means that the 
distance of the pnictogen or chalcogcn to the other 
neighboring Fe layer in turn gets even smaller for the 
11 family. Then this leads to a more efficient role of 
the chalcogen to bridge the two neighboring Fe layers 
to make larger dispersion for the 11 family. On the 
other hand, the distance between the pnictogen and the 
other neighboring Fe layer for the 1111 family is too 
far to make the appreciable three-dimensional dispersion. 

We now summarize the characteristic feature specific 
to the d models with their family dependence. The larger 
hybridization between Fe-3d and p orbitals of pnictogen 
for the 1111 compounds leads to larger extensions of the 
Wannier orbitals over the 11 family for the d model, be- 
cause the Wannier functions of the Fe-3d orbitals more 
strongly mixes with the p orbitals for the 1111 family. 
Still, the nearest- neighbor transfer for the 1111 family 
is not particularly large. Indeed it is even smaller than 
the other families including the 11 family, presumably 
because of large in-plane lattice constant. The largest 
transfer (~ 0.4 eV) appears between the two xy orbitals 
because the orbitals are directed to the direction of the 
neighboring Fe sites, where the amplitude is 342, 315, 
410, and 392 meV for LaFePO, LaFeAsO, FeSe, and 
FeTe, respectively. 

On the other hand, the largest next-nearest transfer 
shows little family dependence. Its value is 0.345 eV for 
LaFeAsO, 0.356 eV for BaFe 2 As 2 and 0.348 eV for FeSe 
and FcTc. These largest transfers are for the hopping be- 
tween the diagonal elements between yz or zx orbitals. 
As a consequence of this trend of the nearest and next- 
nearest neighbor transfers, the 11 family is expected to 
show weaker geometrical frustration effects probed by 
the ratio of the next-nearest-neighbor to the nearest- 
neighbor transfers. 

Although the nesting picture predicts a similar peri- 
odicity of the antiferromagnetic order for all the families 
in contrast to the experimental variation, the contrast in 
the frustration parameter will offer a possibility of ex- 
plaining the variation of the experimental periodicity in 
the magnetic order depending on the compounds from 
the viewpoint of the strong correlation physics. 
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The effective Coulomb repulsion averaged over the Fe- 
3d orbitals for the d model is 2.5, 2.8, 3.2, 4.2 and 3.4 
for LaFeAsO, BaFe2As2, LiFeAs, FeSe and FeTe, respec- 
tively, while the overall bandwidth is around 4.5 eV in 
most cases except for LaFePO (5.1 eV) and LiFeAs (4.9 
cV). This indicates that FeSe is located in much more 
strongly correlated region than the 1111 family, while 
the 122 and 111 families as well as FeTe occupy an in- 
termediate region between the 1111 family and FeSe. In 
fact, FeSe has a substantially larger ratio of U/i ~ 14 as 
compared to 9 for LaFeAsO. An important origin of the 
larger U /i for FeSe is the smaller extension of the Wan- 
nier orbitals ascribed to the smaller hybridization with 
the chalcogen-p orbitals. A smaller Wannier spread yields 
a larger bare Coulomb interaction itself. Another origin 
of the difference is the poor screening in the case of the 
11 family ascribed to the fewer screening channels in the 
absence of La-/ and O-p bands. In fact, these bands are 
sources of efficient screening of the Coulomb interaction 
for the electrons on the iron-3d orbitals in the 1111 fam- 
ily. The distance from the Fe layer to the pnictogen or 
chalcogen elements, namely the h value itself contributes, 
where the larger distance makes the screening weaker for 
the 11 family. 

In other words, an overall trend of the family depen- 
dence of U is explained in the following way: The bare 
Coulomb interaction v is determined predominantly by 
the local environment and the extension of the Wan- 
nier orbital in the real space. On the other hand, the 
screening effect measured from U/v shown in Table IX 
is mainly determined from the band structure in the mo- 
mentum space. In particular the presence of other bands 
close to the Fermi level is important, which efficiently 
screens the interaction of the Wannier orbitals of the tar- 
get bands. As a result, the 1111 family has weaker effec- 
tive Coulomb interaction U as a combined effect of the ef- 
ficient screening by many p bands leading to smaller U/v, 
and the larger extension of the Wannier orbital leading 
to a smaller v. 

Although the overall trend is understood from h, h 
alone is not enough to explain all the details. There 
are some exceptions, and they are ascribed to other fac- 
tors (e.g. number of screening channels, energy differ- 
ence between occupied and unoccupied states). For ex- 
ample, FeSe has larger value of U than FeTe, because the 
chalcogen-p level is deeper (Sec. 3. 4). The U is not much 
different between LaFePO and LaFeAsO. This is because 
difference of the anion (P and As) radius partly cancels 
the difference in h. Larger d bandwidth and higher La-/ 
level weaken screening in LaFePO, which also reduce the 
difference between the two compounds. 

The d model indicates that U strongly depends on 
the orbital for the 1111 family while this anisotropy 
is relatively weaker for the 11 family. In general, the 
x 2 — y 2 orbital has the weakest U because of the largest 
extension of the Wannier orbital extending to the pnic- 
togen site. However, for FeTe, the orbital dependence of 
the extension of the Wannier orbitals is weak because 
of the weak covalency between Fe and Te. Then, U 
for the x 2 — y 2 orbital (~ 3.6 eV) is even larger than 
those for yz and zx orbitals ~ 2.9 eV in contrast to the 



prominently weak U of the x 2 — y 2 orbital for the 1111 
family. 

The larger magnetic ordered moment observed in the 
122 and 11 families may basically be understood from 
the difference in U/i. The larger moment may also 
be ascribed to the lack of the pseudogap structure in 
FeTe, which makes the correlation effect larger, while 
the 1111 family is located in the region characterized by 
"semimetal" with weaker correlation. The exchange in- 
teraction (Hund's rule coupling) J averaged over orbitals 
is 0.39, 0.43, 0.51 and 0.47 eV for LaFeAsO, BaFe 2 As 2 , 
FeSe and FeTe, which has a trend of the orbital depen- 
dence similar to U. 

In the effective models containing chalcogen- or 
pnictogen-p orbitals, namely dp/dpp models, U ranges 
from 4 eV for the 1111 family to 7 eV for FeSe. This fam- 
ily dependence partly comes from the family dependence 
of v, reflecting the family dependence of the extension 
of the Wannier orbitals as seen in Table IX. The orbital 
dependence is much weaker than in the case of the d 
models and more or less isotropic correlation is justified. 
This indicates that the anisotropy of U in the d models 
arises from the individuality of the p bands. Although 
the dp/dpp models are quite complicated to analyze, they 
have several advantages. First, the isotropic correlation 
is convenient when we combine model calculations and 
LDA; in order to avoid double counting of the correla- 
tion effect considered in LDA, we have to introduce the 
so-called counter term to cancel the double counting in 
the model. In the dp/dpp models, since the interaction 
parameters are almost isotropic (orbital independent), 
the counter term of the double counting is very simple, 
i.e., orbital independent. Second, the dp/dpp models may 
describe the antiferromagnetic state more appropriately 
especially for FeTe. There, while the Fe-3c? bands are not 
so isolated, the exchange splitting is expected to be ap- 
preciable since the magnetic moment is 2.2 /lib- Another 
possible advantage of the dp and dpp model is, we can 
explicitly study the polarization effect of the anion-p or- 
bitals, which has been proposed to play an important 
role in the pairing mechanism. 83,84 ) 

Now we have elucidated that this systematic variation 
of the parameters of the effective low-energy models from 
the 1111 to 11 families revealed in the present study is in 
many aspects originated and understood from the varia- 
tion of h, resulting in the change from a strong covalency 
to a strong ionicity in the order from LaFePO, LaFeAsO, 
LiFeAs, BaFe2As2, FeSe to FeTe. The larger correlation 
parameter U/i evolves basically in this order, while the 
frustration parameter t' /t diminishes more or less in the 
same order and the three-dimensionality also evolves in 
this order, all of which are understood from this inter- 
play of the covalency /ionicity. The effect of covalency has 
already been discussed in the comparison of the band 
structure between LaFePO and LaFeAsO by Vildosola 
et al.. 6S > The present work has revealed that this key 
parameter h controls in a more profound fashion the pa- 
rameters of the effective low-energy models including the 
correlation effects and the geometrical frustration. 

Although the correlation effect was not clearly visi- 
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ble in the core level spectroscopy of LaFcAsOi_ 3; F x , 45 - ) 
this can be different for the 11 family in this respect 
of the difference in correlation amplitude. Recent soft 
X-ray photoemission measurement 49,50 ) has elucidated 
the valence band spectra of FeSei-z with x ~ 0.08. In 
particular, Yoshida et at has observed, by the 140 eV 
photon source, a clear peak at the shoulder of the Fe-3ti 
band around 2 eV below the Fermi energy with a dip 
around 1 eV between the shoulder-like peak at ~ 2 eV 
and a rather sharp coherent structure around the Fermi 
level (see Fig. 1 of ref. 49). Although the authors did 
not analyze in detail, this structure does not appear to 
be fully consistent with our density of states shown in 
Fig. 2. Instead, this structure primarily originated from 
Fe-3d bands is reminiscent of the splitting of the coherent 
and incoherent part arising from the electron correlation 
effect. 85 ^ It would be intriguing to perform experiments 
for FeSe and FeTe to see whether the coherent-incoherent 
splitting as well as the satellite structure specific to the 
electron correlation effect exists or not in more detail. In 
addition, although FeSe has the largest U/i in the se- 
ries of the compounds, it is known to be paramagnetic 
metal with the superconductivity at low temperatures, 
while FeTe shows antiferromagnetic order with a large 
ordered moment, naively suggesting a larger correlation 
effect in FeTe rather than in FeSe in contradiction with 
our obtained value of U. It might be related to the fact 
that it is difficult to synthesize the purely stoichiometric 
compound of FeSe. It is important to control the sto- 
ichiomctry experimentally for the purpose of clarifying 
the correlation effect expected in the 11 family. 

In the core level spectroscopy of X-ray, it has been 
claimed that the onsite interaction U is smaller than 2 
e V, 46 ) in apparent contradiction with the present model 
parameter values of the 1111 and 122 family, in the range 
of 2-3 eV. We note that the "U value" speculated from 
the comparison between the photoemission data and the 
small cluster diagonalization without the screening by 
itinerant electrons should be taken with care. In fact, the 
random phase approximation taking into account the full 
polarization including iron-3d electrons generates much 
smaller screened interaction U ~ 0.9 eV, roughly one 
third of the present cRPA estimates. Since the cluster 
diagonalization ignores the metallic screening, the U U 
value" consistent with the photoemission results should 
be similar to this full RPA value, while this fully screened 
"U value" turned out to be a substantial underestimate 
of the model parameter. Therefore, the obtained model 
parameters and the overall photoemission data do not 
contradict each other, but are rather consistent. At the 
same time, it does not mean that the correlation effect 
is small, because the orbital fluctuations and orbital re- 
structuring may be caused by the energy scale of the 
present U obtained by the cRPA for the effective low- 
energy model. If the present interaction by the cRPA 
is comparable or larger than the bandwidth of the tar- 
get bands, we expect substantial correlation effects. This 
is particularly true for the multi-orbital systems with 
screened interactions strongly dependent on the orbitals. 
In iron-based superconductors, this criterion for the sub- 
stantially strong correlation appears to apply. 



The material dependence systematically revealed in 
this study will serve in clarifying the position of this se- 
ries of compounds in the parameter space between the 
weakly correlated systems and the strongly correlated 
region governed by the Mott physics. Although physical 
properties of LaFePO might be understood rather well as 
a weakly correlated system with the nesting scenario, the 
weak correlation picture becomes more and more prob- 
lematic in the other compounds. In the 11 family, phys- 
ical properties can be understood only by considering 
much wider energy range of the excitations away from the 
Fermi level in the order of the bandwidth, if one wishes 
to start from the perturbative picture. The local picture 
of the strong correlation becomes more applicable to the 
11 family. It is remarkable that the ordered magnetic 
moment roughly increases when the effective Coulomb 
interaction increases. The mechanism of the supercon- 
ductivity can also be pursued in this circumstance of 
the correlation effect, where the superconducting critical 
temperature T c appears to be optimized in this interme- 
diate region, namely the region of a strong crossover or a 
possible quantum critical point. An important novel as- 
pect is that the orbital degeneracy and fluctuations of the 
five Fe-3c£ orbitals may play a crucial role in physics of the 
iron-based superconductors in this intermediate region. 
The yz / zx and x 2 — y 2 occupations are largely fluctuat- 
ing in this respect and shows a strong crossover from the 
semimetals to the strongly correlated " half-filled" sys- 
tem with an orbital-selective crossover. Since this region 
is hard to approach neither from the weak coupling nor 
strong coupling limits, these may be studied in the next 
step of the three-stage scheme, where we will continue 
to numerically solve the low-energy models by reliable 
low-energy solvers. 

There is an important caution when one performs 
model calculations by using the parameters determined 
in the present study: The present parameters are derived 
for the real three-dimensional (3D) system and therefore 
the derived parameters should be regarded as inputs for 
the 3D lattice-Fermion model. Conventionally, however, 
iron-based materials have been frequently studied as two- 
dimensional (2D) layer models. If one wishes to extract 
parameters for a purely 2D effective model, one has to 
introduce an additional treatment in the constrained cal- 
culation; when a layer being the target of the 2D model 
is specified and distinguished from other layers, screen- 
ing from the other layers should be included in the po- 
larization calculations, with excluding only the screen- 
ing by the polarizations within the target layer itself. In 
the iron compounds, each layer is metallic, so this inter- 
layer screening effect resulting from the treatment above 
could give nonncgligible corrections of the frequency- 
dependent dynamical screening to the presented param- 
eters for the 3D system. Estimates of this correction are 
important in connecting the ab initio calculations with 
the 2D model analysis in a realistic way, which will be 
presented elsewhere in more detail. 

In the present formalism, the effective low-energy mod- 
els of the d or dp/dpp types have been derived and pro- 
posed by eliminating the degrees of freedom originating 
from the other bands. In this procedure, we have em- 
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ployed LDA/GGA for obtaining the global band struc- 
ture and the cRPA for calculating the screened inter- 
action. In principle, the low-energy electronic structure 
should be determined in a self-consistent way by con- 
sidering the eliminated degrees of freedom again. This 
possible feedback effect has been ignored in the present 
treatment. This is justified under the circumstance that 
the bands of the eliminated degrees of freedom are well 
separated from the target bands near the Fermi level. 
The present iron superconductors appear to satisfy this 
condition rather well, while it is possible that the feed- 
back gives a small but finite quantitative correction. Es- 
timates of this correction is left for future studies. 
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Note added in proof 

Very recently, N. Qureshi et at (ArXiv: 1002.4326) have 
reported that the magnetic moment of LaFeAsO is 0.63 
[lb , which is significantly larger than the values (0.36 
[Ib) previously reported in Ref. 10. This is, however, still 
much smaller than the LSDA estimates. 



